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SUMMARY 


Active magnetic radial bearings are constructed with a combination of permanent magnets to 
provide bias forces and electromagnets to generate control forces for the reduction of cost and the 
operating energy consumption. Ring-shaped permanent magnets with axial magnetization are 
attached to a shaft and share their magnet stators with the electromagnets. The magnet cores are 
made of solid iron for simplicity. A simplified magnetic circuit of the combined magnet system is 
analyzed with linear circuit theory by approximating the characteristics of permanent magnets with 
a linear relation. A linearized dynamical model of the control force is presented with the first-order 
approximation of the effects of eddy currents. Frequency responses of the rotor motion to 
disturbance inputs and the motion for impulsive forces are tested in the non-rotating state. The 
frequency responses are compared with numerical results. The decay of rotor speed due to 
magnetic braking is examined. The experimental results and the presented linearized model are 
similar to those of the all-electromagnetic design. 


INTRODUCTION 


Active magnetic bearings that consist of permanent magnets giving bias flux and 
electromagnets generating control flux have been considered for the reduction of the high cost and 
the running energy consumption in the all-electromagnetic design. Reference 1 reviewed works 
related to this subject, but only a few reports are available to the authors. A key to this subject is 
the construction of magnetic circuits of the combined magnet system. Allaire et al. 1 discussed the 
design, construction and testing of a set that have permanent magnets in the stators. They also 
presented an analysis of the magnetic circuits. Their construction may be advantageous in use at 
high speeds. In addition, it is possible to combine a radial bearing with a thrust bearing. 

However, because the permanent-magnetic bias-flux is not canceled out by the electromagnetic 
flux, the control force may be restricted. 

This paper considers another construction in which the permanent magnets are attached to a 
shaft and their magnet stators are all shared with the electromagnets. This construction has a 
feature that it is possible to cancel out the bias flux with the electromagnetic flux. This point is a 
property similar to the all-electromagnetic design. However, since the composition of the rotor 
becomes complicated, there may be some problems to the rigidity of the shaft and to the elaborate 
machining of the rotor. This magnetic composition seems similar to an all electromagnetic case 
consisting of two independent magnet coils that are wound around the same poles and are assigned 
the separate functions of bias and control. 
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The magnet cores of the experimental setup are made of solid iron for simplicity of 
manufacturing. A simplified magnetic circuit of the combined magnet system is analyzed with the 
linear circuit theory by approximating the characteristics of permanent magnets with a linear 
relation. A linearized dynamical model of the control force is presented with a first-order 
approximation of the effects of eddy currents. The control systems are designed with simple 
analog PID compensators. Frequency responses of the rotor motion to disturbance inputs and the 
motion for impulsive force are tested in the non-rotating state. The decay of rotor speed due to 
magnetic brake is examined. The frequency responses are compared with numerical results. 


PERMANENT-MAGNET BIASED RADIAL MAGNETIC BEARINGS 


Figure 1 shows the mechanical part of the experimental setup of five-axis-control magnetic 
bearing system supporting a symmetric rotor in the horizontal direction. The radial bearings are 
composed of permanent magnets (ferrite) for biasing and electromagnets for control. The 
permanent magnet rings with axial magnetization are attached to a shaft of aluminum and are held 
between two iron rings each. The electromagnets are constructed with a pair of magnet coils that 
are wound around the stators positioned on the radially opposing sides: die pair of magnet coils 
are connected in series and are driven by a single power amplifier. The magnet poles in the vertical 
direction are one-fourth wider than those in the horizontal direction to supply larger bias flux for 
the suspension of the dead weight of the rotor. Four displacement sensors of the eddy current type 
are used to measure the radial motion of the rotor. The magnet cores are made of solid iron for 
simplicity of manufacturing. The rotor mass is 0.90kg, the airgap length is about 0.5 mm and the 
diameter of the rotor core is 40mm. The primary data are given later in Table 1. 

A thrust bearing is positioned in the center of the setup using an electromagnet in place of the 
permanent magnet, for simplicity. Two magnet coils are wound around the inner core of the 
stator, one is for biasing and the other for control. The two opposing magnet coils on the left and 
right sides of the disc (the rotor of the thrust bearing) are connected in series to give similar 
magnetic circuits to the radial bearing. The rotor rotation is provided by compressed air blown 
onto teeth on the edge of the disc. 



Figure 1. Construction of permanent-magnet biased magnetic bearing. 
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Figure 2 gives the magnetic flux paths of the radial bearing, neglecting leakage flux. The 
permanent magnetic flux, providing bias flux, flows radially in the iron-core ring and goes into the 
four poles of the stators through the airgaps, and passes through the stators along the axial 
direction to return to the rotor via the other airgaps. The electromagnetic flux generating from the 
upper and lower magnet coils passes down through the rotor along the radial direction on one side, 
and passes up through the other side. Thus, if the total flux increases in the lower side, then it 
decreases in the upper side. The difference of the flux produce a control force. With this 
construction, we can cancel out the bias flux with the electromagnetic flux. This is a point similar 
to the all-electromagnetic case. The configuration of the rotor, however, may be disadvantageous 
in use at high speeds, and there may be some problems in manufacturing the rotor. 



Permanent magnetic flux ► 

Electromagnetic flux ► 

Figure 2. Magnetic flux paths of radial bearing. 
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Figure 3. Simplified magnetic circuit. 
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ANALYSIS OF MAGNETIC SYSTEM 


Magnetic Circuit of Radial Bearing 


The characteristics of permanent magnets are considered based on the demagnetization curve. 
This curve may be quite linear for permanent magnets useful in practice such as ferrite magnets and 
rare-earth material magnets. Here, we approximate the demagnetization curve with a linear 
relation, using imaginary magnetomotive force and magnetic internal resistance. Then, as a model 
of the magnetic system of the radial bearing we consider the magnetic circuits in Fig. 3, where flux 
leakage is neglected. In the figure, F ^ is magnetomotive force, is magnetic resistance and 
<X>j is magnetic flux associated with a stator numbered j in Fig. 2 and with the permanent magnet 
with B. The imaginary magnetomotive force and the magnetic internal resistance R n ^ may 
be described as 2 

B l 

F mB = , RmB = , ( 1 ) 

Rr M- V ' 

where B r is the residual flux density, [i r the recoil permeability, l m the thickness and A m is the 
area. For this circuit we have the following relations. 

F m\ + Fm2 ~ ^ml^l + F m2^ 2 

F m \ ~ F m 3 ~ ~ ^13^3 

F m \ + F m 4 = R ml O, + R m 4 0 4 ^ 2 ) 

F ml + F mfl ~ Rmft\ * 

% = <*>2 - ^3 + ^4 + 

The magnetic resistances are estimated by 


R. 
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2/ 
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where / mj is an equivalent airgap length considering the magnetic path in the magnet core, m is 
the permeability of air, and A- is the area of a magnet pole. Then, from eq. (2) we obtain the 


relations 
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In the following, we take the center of the bearing as the nominal position of the rotor, and in 
place of the flux we use variables defined by 


Qj 


3nlS.0 

N >’ 


j “ 1,2; Q, 


R* 
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N 


y-3,4 


(5) 
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where /? mj0 are the nominal magnetic resistances, and N is the coil turns that are all the same 
here. The variables have the same dynamical characteristics as the flux and the unit of current, and 
give coil currents if the eddy-current effects are negligible. The variable was introduced first to 
simplify the description of electromagnet systems, and effectively used in modeling the systems 
with eddy-current effects. 4 Then, we rewrite eq. (4) as follows: 

c ml 0, = (l - + « 2 /m2 " U 4 + Qb 

C m202 ~ n \fm\ + (l ~ °2 )fm2 + /m ~Qb (g) 

C m303 ~ ~f\2 + (l ~ )fml + a + Qb 
c m 4 0 4 - fn +a ifm3 + (l - «J La - Qb 

WhCrC ('mB Q/l fn ^ /m ^ ^ /? ) fmB 


C ”V „ ) fmj > J 1*2,3, 4, B 

n mj 0 iV 

fn ~ n \fm\ ~ a 2 fm 2 ’ Jm = a ifml ~ a 4fm4 

We refer the initial steady-state values with subscript 0. In the steady state without initial load in 
the horizontal direction, we have 


Qio ~ Ao + 0/?> 020 = Ao Qb * 

030 = - 04O = Qb > 

mlO = fm2D ~ Ao* f m 20 ~ fm4Q ~ Ao = ® 


(7) 


In the following , we assume that 

(1) The displacement of the rotor is sufficiently small. 

(2) The incremental magnetic flux is sufficiently small. 
We consider the increments 


Q, = 0/ “0, 0. A fry = /„ ~ fry 0 * 

A a. «a y -a yo (8) 


We write the magnetic resistance as 


^ml - ^mlo(l + ^ml)» ^m2 ^mloO ^ml)> 

^m3 = ^b3o( 1 + ^m3)> ^m4 = ^m.3o(l _ Al) 


(9) 


where is the relative variation of the magnetic resistance. Then, the parallel magnetic 
resistance R T defined under eq. (4), is expressed as 

R t * Bnjl- 2(« ]0 C + cc^kj)] * B ro 

where a ]0 and ot^ are the nominal values. This implies that the parallel resistance is nearly 
constant under the assumption (1). This gives the increments of a, as follows: 


-Am, =■ tia 2 = n Tw k ml , - A« 3 = A a 4 = ajc m2 


(ID 
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Substituting eq. (1 1) into the incremental equations of eq. (6), and neglecting the second-order 
terms of the increments, we obtain the following relations. 


where 


“ (l **lo)4/ml 0f 3o4fm34> 

” + ~ a io)4/m2 + **3()4/m34 

(12) 

“ (l “ ^30 )4/m3 + **3<)4/m4 “ 

- OIjQ + (l “ ^30 )4^m4 * **10 4/|2 

(13) 

" 2« 10 / 10 (l 2 c* 10 )/jq +■ 

) + 2flq 0 /|o = (l 2020)^10 6b 

(14) 

- 2/ 10 * mJ , 4/ ml2 = Af mX - Af ml , 

~L 4 

(15) 


Dynamical Linearized Equation of Incremental Flux 


We discuss the vertical direction whose special case gives the results of the horizontal 
direction. We apply a model presented in Reference 3, considering the first-order effect of eddy 
currents. This model describes the electric and magnetic coupling system composed of the two 
stator-coils in series as 


V^R n l^N(<P^0 2 ) (16) 

(/y e + Lj e + N4>). = 0, (RJ L + =0, j =1,2 (17) 

^ i-1,2 (18) 


where V,, R n and /, are the input voltage, resistance and current of the magnet coil; i ej and i u 
are the eddy currents in the magnet cores, R eX (= R t2 ) and R u (= R L2 ) are the resistances, and L el 
(= L e2 ) is the inductance of the equivalent eddy current loop; are the magnetomotive forces. 
Here, we take the coil current positive when it increases the flux in the stator 1 (in the lower side). 
For a power amplifier with output current control, we suppose that the input voltage is expressed 

as , 

v x = rib E l - /,) (19) 


where E x is an input voltage to the amplifier, p is the control loop gain, and b ’ is a constant. In 
this case, eq. (16) is replaced by 

pb' E, - [p + R i2 )/, + iv(*, + # 2 ) ( 20 > 

Using the variables defined by eq. (5) simplifies eqs. (17) and (19) as follows: 


/, -6,£ 1 -7„(e l +a) 

(*"e + TJ e ) “ -TM, ijj m ~T ol jQj, j = 1,2 


( 21 ) 

( 22 ) 
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where 
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Eliminating I v i ej and i Lj from eqs. (18), (21) and (22), we obtain the following equations. 


where 


+ (Tj, + E , ) - T n Q, -T„Q,- T„(t.& + ft ), 
tJ m *U- *e,)-t 2 ,q 2 - t„q 2 - +a) 

^21 = + T ol \ T eV T u — (t r + T ol + T(j e \ 


(24) 


The interacting magnetomotive forces in eqs. (12) and (13), A/ ml2 and A/ m34 , which are zero 
if there is no effect of eddy currents, make the analysis quite complicated. We suppose that these 
influences are secondary, and we will make a rough approximation. From the incremental 
relations of eq. (24), we have 

7«l4/ml2 4 /" ti 12 = ~ ^0il^el9l2 ~ > *?12 ~ Q\ ~ * 


T 0 u t T oia + T ( 


Oel 


(25) 


This relation suggests that the interactions are dynamical, i.e. statically zero. On the other hand, 
eqs. (12) and (13) give the relations 


where 


^12 = 012 k m \ + (1 2tt I0 )/y ml2 2 n 3 a^fn 34 > 

9.34 = ^^lofio^mi + (1 “ 2r* 3n )-4/’ m34 - 2« 10 4^ nl2 > ^34 = Q-s ~ q 4 


(26) 


012* = 0, o' <?2o' = 2(l-2« IO )/ IO 

Substituting eq. (26) into eq. (25) and a similar equation of A/^, and taking the lowest derivative 
terms of the variables, we have the approximation 


4/m 12 02 Wm„ (27) 

4/m 34 = 


We eliminate A/ ml and A/j^ from eq. (12) with (27) for A/ m34 and the incremental equations of 
eq. (24). Since the results have higher derivatives of k ml and are so complicated, we neglect the 
higher derivatives than the first. Thus, we obtain the following results. 


+ + +(r, + 0 , 0 * + 0 ,o'^m, +eX 

^2i^2 + T n q 2 + q 2 + (-y, + 03 d’ rj* ml + CV*-, “ +<?, ) 

where 

4 - (2J-,, t T„X , f„ - 2T„ + r m , + + r„, 

^1 " 2/ 10 [(l -2& w )T Rl - 2«, oOfjo 7^4*3 ], 7^,'= 7" ri +«, 0 (7’ oil + 7^,) 


(28) 

(29) 


The time constants 7^, and T u are the same values that are given when we take the two 
opposing stator-coils to be one. This natural fact is easily confirmed with the definition of the time 
constants in eq. (23). Thus, the above results are similar to those of the all-electromagnetic cases. 
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except for the term with Yi- This term appears without eddy current effect. We can see that its 

first term with To,’ results from the effect of the opposite-side stator-coil system with a bias 

current We note that this term is reduced by the second term that is due to the interaction from the 
horizontal direction. 

Similar equations are derived for the horizontal direction, but it may be simpler to take another 
way. Adding the two equations of eq. (13) gives the equation without interaction to the summing 

variable as 

2 <730 + 2Q£ m > = 4^3 + 4f«4 > 2<730 = <73 + <74 (30) 


Then, using relations similar to the incremental form of eq. (24), we have 

7^30 + T„q 30 +<730 + Qb( IKi + 0 “ + e a ) (3 1 ) 

In a similar way, the difference variable q M , with A/ ml2 in eq. (27) we obtain the approximated 

(1 - 2 « 3fl ) 7 ' ot3 r e 3 < 7 34 +[(1 - 2« 30)7^3 + 7 ^j] <734 +<734 

= 4« 10 / 1 o{[T rf -(l-2« 1 o)7' 0 „ ] ]fc ml +* ml } (32) 


Then, we have the expression of the decoupled form 

*73 “ *?30 + *?34 / 2 , <74 = *730 _ *?34 ^2 


(33) 


Similar to the case in the vertical direction, the interaction due to the term A f mn reduces the effect 


The above results suggest that the interactions are primarily due to the bias current combined 
with the airgap variation in the biasing direction, and that this interaction into the horizontal 
direction becomes smaller at higher frequencies. The results become quite simple in the case of 
laminated magnet cores. 

We will obtain the time constants in eqs. (28) and (31), experimentally, because it is difficult to 
estimate theoretically the time constants of the eddy current effects, defined in eq. (23). For the 
fixed airgap of the nominal length (& ml = k ml = 0), we write eq. (28) with the parallel combination 
of two first-order time-lag systems, using the Laplace transforms, as 

Q(s) fo + 1 *1 , l-*i 

bE(s) t 2 s z + 7^+1 T v s + 1 T 2 s + 1 (34) 


where we omitted the subscript 1, and where 

T 2 = T x T 2 , T x ~T x + T 2 , T t = (1 - k x )T x + k x T 2 

To give Yi > the other time-constants, T ou and T M , are calculated from these results and the 
theoretical value of T R1 with the relations in eqs. (29). 
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LINEARIZED EQUATIONS OF ROTOR MOTION 


Linearized Equation of Control Force 

The magnetic force acting between two pole faces with the magnetic flux <t> may be expressed 

as F = ( AJ 2p 0 )( O / A ) 2 , where A is the area of the pole face and A h is its projected area on 
the plane perpendicular to the sense of the force. Applying this to the two poles with the same face 
area and using the variables of eq. (5), we have the resultant magnetic force acting on the rotor on 
the opposing sides as 


Fn - c F ,(Q, 2 -G 2 ! ). 7^ 

(36) 

Then, the net force is given by 


= fal ~ ^120 ~ f\ +C Fi ) 

(37) 

with the first-incremental force 


fa ~ (^fi/2)(<7i 

(38) 

where 


*F| - 4c fl 0, - 4 10 , c o ,--^(>0) 

(39) 


When the rotor approaches stator 1 with displacement z, the relative variation of the magnetic 
resistance, k mi , may be approximated under the assumption (1) as 


k 


ml 



(40) 


Using eqs. (28) with this relation, and replacing ( q l + c Ql q 2 )/ 2 by a new variable q x , we write eq. 
(31) with the dynamical form as follows: 


where 


fa T 2l q t + T n q x +q x -a n z-a 0l z = f\b x {T eX e x +e,) 

a n = +£, ). a m = ^-(l + c G1 c ei '), 


10 


21. 

(j = z —(\ + c ) c ' _ ^ 2 ° 

f » 7 \ 1 + c ei/> c q \ ~ 777 » 

*£10 


T'=— C ° 1 ^ 1 Ssl2(t 1 ') _ 2a, 


1 + Cq x Cqi 


10 


‘(X T 

1/1 * n i 


1-2 a 30 ' (tU3 ’ 

1 


(41) 


(42) 


For the horizontal force, the result is directly obtained by application of eq. (31) to eq. (38) with 
c Q1 =1 as follows: _ q 

fa ~ + faifa + Qi ~ 7 ^“( + y) - b 3 {T e3 e 3 +■ e 3 ) 

l m30 (43) 


where y is the displacement toward stator 3. 
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The above results are almost all the same as those of the all-electromagnetic case, 4 if we regard 

the equivalent bias flux as an equivalent bias cun-enL The term with TLS is inherent in the 

magnet system consisting of two opposing stator-coils with bias current. We see through the 
derivation that the interaction appears separately on each of the two opposing sides, but that their 
resultant effect becomes smaller, and is canceled out in the horizontal direction. In this way, our 
analysis may seem to be a waste of labor and space. The analysis, however, gives a deeper 
understanding of the phenomena, and shows the complexity of the magnetic system. 


Linearized Equations of Rotor Motion 


The radial position of the rotor is controlled by the resultant forces and resultant moments of 
the incremental forces in the two bearings. We number the bearings 1 and 2 on the left and right 
sides respectively, and use these numbers in the subscript of the variables if necessary , tor 
simplicity We take the sense of the rotor displacement to be positive when going down, as in the 
above and the sense of the conical motion to be positive when going down more in bearing 1 than 
in bearing 2 We rewrite the conical motions with the equivalent displacements at the bearing, tor 
convenience, and neglect the interactions due to the gyroscopic effects. Applying eq. (41) to the 
two radial bearings having the same characteristics and supporting a symmetric rotor, we obtain the 
linearized equations that are similar in form each other, for the rigid-rotor motion m the vertical 
direction, as follows: 

Translatory motion: mi' - *,A + 4, 


Conical motion: 


T 2l q z + f u q z +q z - a n z-a 0l z = PM T <A + 0 

m <Ac “MzC +</ zO 
^21 tfzT. + + < 7 z c ”^01 


(44) 


(45) 


where d is disturbance force, and m c and d lC are the equivalent mass and force of the conical 

5 Z r™ : ^ a ^ tr\ thp actual rnntrnl inrmt& as 


wnere u z is uisiuiuaiii'C iuicc, auu ui c auu ^ zC 

motion. 5 The control variables e t and e zC are related to the actual control inputs as 


2e z ” e n + e 2v le^ = e n e 2l 


(46) 


If the control variables are given, then the inputs, e n on the left side and e 21 on the right side, 
are determined by 


e. + e 


zO 


e 2 i Ctf 


(47) 


Similar equations of the horizontal direction are obtained with c Q1 = c Q1 ’ = 1 in the above results. 


EXPERIMENTAL AND NUMERICAL RESULTS 


The primary data of the experimental setup are shown in Table 1. We mounted the setup by 
pressing on to a block of modelling clay. Experimental data are obtained in the non-rotating state 
except for the decay of rotation. We will show the results in the vertical direction; similar results 
are obtained in the horizontal direction. 
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Table 1 Data of Setup 


Rotor 



Mass (PM contained) 

m = 0.90 

kg 

Diameter 

D = 40.0xl0' 3 

m 

Moment of inertia of conical motion 

2.31x10 3 

kgm 2 

Polar moment of inertia 

J^ = 3.51xl0' 4 

kgm 2 

Distance between bearing centers 

lj = 111 x 10' 3 

m 

Distance between two gap-sensors 

l 2 = 174xl0' 3 

m 

Stator (Radial) 



Magnet coil: Turns 

II 

to 

8 

Turns 

Resistance 

R = 0.62 

n 

Area of pole face (Vertical) 

A y = 1.60x1 O’ 4 

rv » 2 

m 

Area of pole face (Horizontal) 

A„= 1.28x1 O' 4 

m 


/„ = 0.5x10 3 

m 

Thrust bearing 



Bias coil : Turns 

* 

ll 

8 

Turns 

Bias current 

A, = i.o 

A 

Control coil: Turns 

* 

II 

O 

O 

Turns 


Magnetic Characteristics 


It is difficult to theoretically estimate the time constants of our model. We examine the 
dynamical generation of the incremental flux by obtaining its frequency response to an input into 
the power amplifier. The generating voltage in a search coil gives the frequency response of the 

rate of the flux. We have the response of the flux by numerically multiplying l//co to this result 

represented in complex number, where j is the complex factor and to is an angular frequency. 
The result measured with a 4-tum search coil is shown in Fig. 4 by the solid lines in the vertical 
direction, for the fixed airgap length of 0.5mm and the non-biasing input of the amplitude 
corresponding to the coil current of 0.38A st atica lly. The phase in low frequencies is considered to 
be incorrect (the authors guess this is due to FFT analysis in the analyzer). The response has a 
characteristic that the phase lag in higher frequencies is small compared with the decay of the gain. 
We may understand this result with a first-order time-lag model in such a way that the increase in 
the magnetic resistance, due to the eddy current effects, decreases the incremental flux, but reduces 
the time constant that is inversely proportional to the resistance: the ratio of the inductance to the 
resistance. 

To approximate the characteristics, we apply the model of eq. (34), and try curve fitting by 
selecting the two time constants and the weighting. The fitness depends on the range because of 
the low-order model. A result is shown in Fig. 4 by the broken lines; with this curve, we have T x 
= 0.17ms, T 2 = 1.88ms and Jfc, = 0.40. The static gain of -84.4dB is given by the theoretical value 
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calculated with b { = 0.36 and / m0 = 1.05 x (airgap length) without flux leakage nor fringing 
effect The lines zigzag through the experimental results, and may be unsatisfactory: the gain is 
larger in frequencies of about 15 to 100Hz, and the phase lag is smaller as a whole; however it 
seems difficult to have much more improvement. A similar result was obtained in the all- 
electromagnetic case. 6 



Figure 4. Frequency characteristics of incremental flux. 


Control System 


The control systems are constructed as in Fig. 5, for each of the vertical and horizontal 
motions, neglecting the interactions of the conical motions due to gyroscopic effects. This 
construction decouples the control of the translatory and conical motions into two single variable 
systems within the linearity. The input variables z x and z 2 are displacements of the rotor on the 
left and right sides. The outputs e x and e 2 are control inputs given to the electromagnet systems 
on the left and right sides. The variables z and z c are displacements of the translatory and conical 
motions, u and u c are control variables, and C( s) and C c (s) are transfer functions of the 
compensators. The inputs u 0 and are used for the measurement of frequency responses of 
the control system. We used analog PID compensation. 

Gains of the PID controllers were adjusted experimentally to give desirable characteristics of the 
control system. The selected gains gave the phase margin of about 35 deg. at 80~90Hz, and the 
gain margin of 15-18 dB at around 400Hz in the open-loop characteristics. The control inputs e x 
and e 2 , including the bias input, are limited to a value corresponding to the static coil current of 
about 1.5 A. 


732 




Uo 



Compensator Limitter Power amp. 


Figure 5. Control system. 


Frequency Responses of Rotor Motion 

Frequency responses of the rotor motion are measured for a sinusoidal input of amplitude 
0.5V, corresponding to the coil current of about 0.1 9 A statically. Figure 6 shows the responses of 
the vertical 



Frequency [Hz] 


Figure 6. Frequency responses to translatory input. 
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Figure 7. Frequency responses to conical input. 

motion to the disturbance input of the translatory system. The phase is given only for the 
translatory motion. The maximum displacement in frequency about 20Hz is about 0.13mm, about 
one-fourth of the airgap length, (the output gains of the displacement sensor are 5.0 x 10 3 V/m for 
the translation, and 7.8 x l(r V/m for the conical motion). The induced conical motion is small as 
expected, but may have some information in higher frequencies. The responses to the input into 
the conical system are given in Fig. 7. The results are similar to the above results with a small 
translatory motion induced. The maximum displacement is about 0.13mm. 

It was difficult to measure the permanent magnetic flux density in the airgap 0.5mm, therefore 
we guessed roughly 0.2T from the other measurement data. This flux density corresponds to an 
equivalent coil current of 0.84 A ( Q h = 0.84) and estimates the forces of 5.7N and 4.6N in the 
vertical and horizontal directions, respectively. In this case, the bias coil current in the vertical 
direction is calculated as 0. 19 A (/ 10 = -0.19). The data necessary for the numerical analysis and 
the time constants obtained from the experiment, excluding data given in Table 1, are summarized 
in the followings. 

b x « 0.36 A/V, p * 140 V/A 

« I0 =0.25, = 0.20, a w = 0.11, R ml0 =5.0 xlO 6 A/Wb 

T R1 = 0.057 ms, T aLl =0.26 ms, T 0r} =0.82 ms, T el =0.85 ms 
a 01 = 2.11 x 10 3 A/m, o„ = 1.83 As/m, c Q1 = 1.59, c 01 '= 1.26 

We use T 0Lel in place of T 0Le3 in eq. (42) because we have no data for it at present; however the 
authors suppose that this approximation is unimportant in this case. 

The broken lines in Figs. 6 and 7 give the numerical results based on the linearized model. The 
results fit fairly well with the experimental results in lower frequencies, but poorly at frequencies 
higher than 150Hz. A further examination is required to explain this inconsistency. We guess at 
present that the primary cause is in the stiffness of the stators and in the supporting of the 
experimental setup. 
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Motion for Impact 

Figure 8 gives the translatory motion and the actual control input for vertical impulsive force 
added at the center of the rotor for the maximum displacement of about 70% of the airgap. The 
shapes of the transient responses were similar for smaller motions, which is a result of the lower 

control gains. 



Figure 8. Translatory motion for vertical impact. 


Decay of Rotation due to Magnetic Braking 

We had a maximum rotating speed of a little more than 7,000 rpm, which was restricted by 
tolerable whirling of the rotor. The decay of the rotor speed starting from 7,000 rpm is observed 
as in Fig. 9. The rotor stopped completely after about 55s. A line connecting the points is 
approximated by the equation ci> + c tw = -r b , c - 0.046 rad/s, r b - 72 rad/s 2 . 



& 



0 10 20 30 40 50 60 


Time [s] 

Figure 9. Decay of rotation. 
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CONCLUSIONS 


Active magnetic radial bearings were considered with the combination of permanent magnets to 
give bias force and electromagnets to generate control force. The ring-shaped permanent magnets 
with axial magnetization are attached to a shaft and share all the magnetic poles with the 
electromagnets, so that the control flux may cancel out the bias flux. All the magnet cores are made 
of solid iron for simplicity of manufacturing. This, however, caused much complicated analysis. 
The simplified magnetic circuits were analyzed and the linearized dynamical model of control force 
was presented with the first-order effect of eddy currents. Frequency responses of the rotor 
motion to disturbance inputs and the motion for impulsive force are tested in the non-rotating state. 
The frequency responses are compared with the numerical results in some disagreement. The decay 
of rotor speed due to magnetic brake was examined. The experimental results and the presented 
linearized model are similar to those of the all-electromagnetic design. 

We considered larger, but the same-size magnet poles in the vertical direction supporting the 
deadweight of the rotor, to reduce a decrease of the resultant bias force on the lower side. It may 
be reasonable to design different-size poles so that the permanent magnet may supply different bias 
fluxes to support a part or all of the initial load. In this case, however, the analysis of the magnetic 
circuits becomes much more complicated. Our approach may be a step towards this case. 

The authors thank Mr. S. Fujino for help in making the experimental setup. 
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